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Abstract. 

This paper is firstly intended to review shortly the most recent developments and ideas resulting from the necessity of hav- 
ing a scale at which gravitation to unify the other fundamental forces. With the declared intention of predicting an in situ 
distinct possibility for Kaluza-Klein gravitons forma tion in the supernova shell we are using state of the art simulations 
lLanger et al. 19991: luclboula & Owocki 2002t lud-Doula & Owocki 2003k lud-Doula 20041) for massive stellar winds to infer 
that the supernova shock which will hit such winds will meet a non-isotropic and non-homogeneous matter distribution with 
a very distinct geometry. By linking this to the mechanism of particle shock acceleration at cosmic ray energies (Fermi accel- 
eration) and the related spallation in the wind shell, the result is the creation of proper conditions for Kaluza-Klein gravitons 
formation in the supernova shell from neutrino secondary particles (cosmic ray spallation products) interacting with ~ 10 18 eV 
cosmic rays. 
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1. Introduction: Brane-World Scenarios 

The main motivation for considering the spacetime as multi- 
dimensional, by introducing extra-dimensions, comes from the 
need to provide an unified theory in which two separate scales, 
the electroweak M w ~ 100 GeV and the Planck scale M P ~ 
10 19 GeV can coexist (hierarchy problem). The string theory 
and M-theory are able to incorporate gravity in a reliable man- 
ner. Recently a new multi-dimensional the ory splitted from 
string theory: brane theory. From Nords trom ( Nordstrom 1914) 
and Kaluza-Klein JKaluza 192 it iKlein 1926al iKlein 1926bl) 
initial five-dimensional scenarios, with one circular, compacti- 
fied dimension and the constraint that the field must not depend 
on the extra-dimen sion, where developed three main classes o f 
extra-dimensions (Hewett & Spiropulu 2002; Rub akov 20011) : 



1. Large extra-dimensio ns 

Arkani-Hamed et al. 19981: 



Dick 2001) 



(Antoniadis et al. 1998; 
lArkani-Hamed et al. 1999 : 



The Standard Model gauge and matter fields are confined 
to a three-dimensional brane that exists within a higher 
dimensional bulk. Only gravity can propagate into the 
compactified extra spatial dim ensions. 
2. Warped extra-dimensions jRandall & Sundrum 1999at 
iRandall & Sundrum 1999bl) 

The hierarchy between the Planck and the electroweak 
scales is generated by a large curvature of the extra- 
dimensions. The simplest framework in this scenario 
is is that of a five-dimensional Anti-de-Sitter (AdSs) 
space, which is a space of constant negative curvature. 



In the compactified, finite extra-dimension just gravity 
propagates. 

3. TeV~ l -Sized extra-dimensions jAntoniadis 1990l) 

By themselves, they do not allow for a reformulation of 
the hierarchy problem, but they may be incorporated into 
a larger structure in which this problem is solved. In these 
scenarios, the Standard Model fields are phenomenologi- 
cally allowed to propagate in the bulk. This presents a wide 
variety of choices for model building: 

- all, or only some, of the Standard Model gauge fields 
exist in the bulk; 

- the Higgs field may lie on the brane or in the bulk; 

- the Standard Model fermions may be confined to the 
brane or to specific locales in the extra-dimension. 
The phenomenological consequences of this scenario 
strongly depend on the location of the fermion fields. 



2. Brane Formation by Particle Collision 

In principle, higher-dimensional objects can be formed also in 
particle collisions (p-branes). The standard picture is that when 
the colliding partons have a center-of-mass energy above some 
thresholds of order of the Planck mass (~ 10 19 GeV) and the 
impact parameter less than the Schwarzschild radius, a Black 
Hole (type-0 brane) is formed and almost at rest in the center- 
of-mass frame. The Black Hole so produced will decay ther- 
mally and thus isotropically to that frame. In a universal extra- 
dimensions scenario, Kaluza-Klein (KK) states of quarks and 
gluons could be produced in colliders. 
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In the large ext ra-dimensions scena ri o the first class of col- 
lider processes JGiudice et al. 1999t iMirabelli et al. 19991 
ICheung & Keung 19991) involves the real emission 
o f KK graviton states in the scattering processes 
( faewett & Spiropulu 2001 : see Fig.[TJ: 



-> y + G {n) otZ + G 



pp^g + G w ; pp^g + G 



in) 



Z^>ff + G (n) 

The produced graviton behaves as if it were massive, non- 
interacting, stable particle and thus appears as a missing energy 
in the detector. 

At hadron colliders, for the interaction qq — > g + G^ is 
espected an emission process which to results in a monojet plus 
missing transverse energy signature. 





Fig. 1. Some interactions which can give rise to emission of 
real KK gravitons. 

Deviations in cross sections and asymmetries in Standard 
Model processes as (see Fig.|2}: 

<? + e~ —* ff or yy or ZZ or WW 

PP -> // or yy 
ep — » e + jet , 

or new production processes which are not present at 
three-level in the Standard Model as, for example gg — > 
ff, might be an evidence for a virtual graviton exchange 
(Giu dice et al. 1999t iHewett 1999 ) but also putting limits on 
the extra-dimension size. 




f V 



KK 



rq V TV 

Fig. 2. Virtual KK graviton exchange possible processes. 

In the warped extra-dimension scenario, the principal 
expected collider signature should be the direct resonant 
producti on of the spin-2 states in the KK grav iton tower 
jDavoudiasl et al. 2000t faewett & Spiro pulu 20021) . 



In the TeV~ l scaled extra-dimensions scenario, the masses 
of the excitation states in the gauge boson KK towers de- 
pend on where the Higgs boson is loc ated. Precision elec- 
troweak data place strong constraints jRizzo & Wells~2 000 ; 
Masip & Poma rol 19991 iMarciano 19~99h on the mass of the 
first gauge KK excitation. 

For a detailed image of the physic al processes in w hich are 
involved the KK states see Han et al. faan et al. 19991) . 

3. Extra Dimensions Astrophysical and 
Cosmological Constraints 

Astrophysical and cosmological considerations impose strict 
constraints on some theories of extra-dimensions. A part of this 
constraints are resumed in tabled 

All the fundamental mass scale limits of multi-dimensional 
gravity from tabled apply only for a situation where all extra- 
dimensions have the same compactification radius. In the gen- 
eral case the bounds could be less restrictive. 

Ia Supernovae , Cosmic Microwave B ackground (CMB) 
and clusters data ( Avelino & Martins 20021) . also can impose 
restrictions on brane world scenarios through the values of 
the cosmol ogical parameters. As we ll, cosmic ray exper- 




iments JAncho rdoaui e \ al. 2' 
Anchordoaui & Goldberg 2' 
are imposing restriction through 
neutrino penetrating showers, 
Standard Model rate, initiated 



iFeng & Shapere 2002} 
lEmparan et al. 2002 ) 
the absence of deeply 
at rates above the 
in the Earth's atmo- 



sphere by brane form ation and decay (vN — » p-brane 
(Kolbe 



; et al. 2003; Langanke & Martinez-Pinedo 2004, 
lLanganke & Martinez-Pinedo 2004bl) . 

4. Size of Extra-Dimensions 



In the toroidal compactification assumption, in the case in 
which there is only one extra-dimension, for large extra- 
dimensions KK massive states there is an infinite-dimensional 
symmetry (the loop algebra on S 1 ) at the Lagrangian 
level dDolan & Duff 1984 lAnlakh & Sahdev 198.4 
ICho & Zoh 19921) . but it is broken by the vacuum config- 
uration. 

Han et al. paper's Introduction JHan et al. 19991) contains 
an important affirmation: "A torus compactification is per- 
haps not realistic, since the bulk fields which we are ignor- 
ing are potential sources of n-dimensional curvature, as are 
the branes themselves. However the torus has the great advan- 
tage of conceptual and calculation simplicity". So, just because 
its "s implicity", many of the authors, with some exceptions 
dKaloper et al. 2000I) . in the field assume that the compactifi- 
cation is toroidal. 

Also, in AD D jArkani-Hamed et al. 1998: 
lArkani-Hamed et al. 19991) multi-dimensional scenario, the 
four-dimensional Planck scale is not a fundamental parameter. 
Rather, the mass scale of multi-dimensional gravity, which we 
will denote by M, is fundamental, as it is this latter sca le that 
enters the full multi-dimensional action (Rubak ov 200 ll) : 



S = - 



1 



16ttG 



(D) 



I 



d D X 



(D) R (D) 



(1) 
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Table 1. Astrophysical and cosmological constraints on the fundamental mass scale of multi-dimensional gravity (M > the table 
value in TeV; d is the number of extra-dimensions) (Hewett & Spiropulu 2002). 
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1) Cosmic SNe 


84 


7 




(Hannestad & Raffelt 20011 


2) yi> Annihilation 


110 


5 




(Hall & Smith 19991 


3) Re-heating 


167 


21 


5 


1.5 (Hannestad 2001) 


4) Neutron Star Halo 


450 


30 




(Hannestad & Raffelt 2001) 


Overclosure of Universe 


6.5/ 






(Hall & Smith 1999) 


Matter Dominated Early Universe 


85 


7 


1.5 


(Fairbairn2001) 


Neutron Star Heat Excess 


1700 


60 


4 


1 (Hannestad & Raffelt 200 1 ) 



where: 



1 



1 



T{D) 



M D-2 M d + 2 



is the fundamental D-dimensional Newton's constant, d — D - 
4 is the number of extra-dimensions, and d D X = d 4 x d d z. 

The long-distance four-dimensional gravity is mediated 
by the graviton zero mode whose wave function is assumed 
as homogeneous over extra-dimensions. Hence, the four- 
dimensional effective action describing long-distance gravity is 
obtained from equation by taking the metric to be indepen- 
dent of extra-coordinates z. The integration over z becomes 
trivial, and the effective four-dimensional action: 



'«// = 



16ttG 



(D) 



J d 4 X y[^R (4) 



(2) 



where Vd ~ R is the volume of extra-dimensions. From eq.|2 
we see that the four-dimensional Planck mass is equal to: 



Mpi =R--M- 



(3) 



Therefore, if R is large then M can be much smaller than 
Mp\. If this scen ario is to solve the hierarchy problem then M < 
10 - 100 TeV (lHannestad & Raffelt 200 It Mannestad 200 ll) . 
Because of this requirement, d — 1 is excluded, M a 100 TeV 
corresponding to R 10 8 cm. 

But, if we will not consider the simplifying remarque that 
the graviton zero mode is homogeneous over extra-dimensions, 
but instead the graviton having an organized and structured in- 
ternal distribution (as we will see further in D.E.U.S. model) 
and, even more, that metric and the Ricci scalar are not inde- 
pendent of extra-dimensions, the integration over z in equation 
will not be so trivial anymore. In consequence, in the four- 
dimensional effective action from equation (f2} will enter sup- 
plementary terms, equation Q becoming unrealistic and un- 
able to describe in a right manner the extra-dimensions radius. 
The result is that the one-dimensional extra-dimension case 
won't be possible to be excluded in the above trivial manner. 

In warped extra-dimensions scenarios (with emphasis on 
the Randall and Sundrum model) the effects of the branes on 
the bulk gravitational metric are taken into account, and are 
expressed in the action by the vacuum energy terms V V i S a,i e and 

* hidden ■ 



S — S gravity S visible S hidden 



(4) 



with: 



' gravity 



f d 4 x f dz yj-g^l-A + 2M 3 R (5) ] 



$ visible ~ ^ d X 'sj ^visible [^visible 
J d 4 X ^~ghidden[Lhidden 



], 



* hidden 



visible 



Vhidden\ 



where A is the cosmological term, and the fifth-dimension bulk 
is located between < z < z c , z c = nr (r is the bulk radius). 
At each boundary plane it is assumed a Z2 reflection symmetry. 
This results in a jump in the extrinsic curvature at these planes, 
yielding two domain branes of equal and opposite tension. One 
ends up with two three-brane rigidly located at the S 1 /Z2 orb- 
ifold fixed points in the z direction while extending in the x^ 1 di- 
rections. The solution to the five-dimensional Einstein's equa- 
tions is required to respect four-dimensional Poincare invari- 
ance in the xf directions, and is found to be: 

ds 1 = e- la(i) ^ v dx* dz v + dz 2 , 

with cr(z) = k I z |. The AdS 5 curvature scale k is the single 
scale with fine-tunes between A, Vvisibie and Vudden, through 
V hidden = -Vvisibie = 2AM\ and A = -2AM\ 2 . 

The relation between the Planck scale in four dimensions 
and the five-dimensional Planck scale M is found by integrating 
out the extra-dimension coordinate in the action: 

= -«-**). 

Unlike the ADD scenario's equation 0, the two Planck 
scales in the Randall-Sundrum scenario are of the same order. 

5. Acceleration and Branes in Supernova Shell 

5.1. Massive Star Winds 

Into the paper ( Popesc u 20041) we see that for massive stars 
with strong winds (3OM < M < 5OM ) the radiative ac- 
celeration plays an important role in the phase space dis- 
persion of different wind particle species before injection in 
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the supernova shock. Taking into account that for a mas- 
sive star the degree to which the wind is influenced by 
the magnetic field can be described by just one dimension- 
less paramete r: the wind magnetic con finement parameter, 77, 
Jud-Doula & Owocki 200*21 lud-Doula & Owocki iool . mean- 
ing the ratio between the magnetic energy density and the ki- 
netic energy density. ud-Doula & Owocki simulations (see Fig. 
13 show that, for a standard O star, with 77* < 1 the magnetic 
field is fully opened by wind outflow, but that for 77* =0.1 it 
can enhance the matter density near the magnetic equator. For 
stronger confinement, 77* > 1, the magnetic field remains closed 
over limited range of latitude and high above the equatorial sur- 
face, but eventually is opened into nearly radial configuration 
at large radii. Within the closed loops, the flow is channeled to- 
ward loop tops into perpendicular shock collisions. Within the 
open field region, the equatorial channeling leads to oblique 
shocks at the magnetic equator. Their assumption of isothermal 
hot-star wind is valid for the dense winds with the characteristic 
wind column densities much bigger than 7 x 10 21 cirT 2 . For less 
dense winds with the columnar densities much smaller than 
7 x 10 21 cm' 2 they used the adiaba tic approxim ation. The sim- 
ulation dud-Doula & Owocki 20021 hid-Doula & Owocki 2003b 
shows that, for 77, < 1 and identical stellar wind parameters, the 
solutions of isothermal and adiabatic models are similar. A dif- 
ferent behavior can be seen for 77, > 1, for which the adiabatic 
models exhibit an greatly extended magnetically confined re- 
gion compared to their isothermal counterparts. These regions 
are filled in with very hot (up to 10 8 K) low density gas that 
is in hydrostatic equilibrium. This is in contrast to isothermal 
confined region where hydrostatic atmospheres are not possible 
due to relatively low sound speed, and stellar gravity pull of the 
compressed, stagnated material within closed loops into an in- 
fall back onto the stellar surface through quite complex flow 
patterns. For the isothermal cases with zero rotation, unlike in 
the adiabatic models, there is a compressed dense, slowly out- 
flowing equatorial disk (see Fig.|3}. Instead, the region above 
the magnetic closed loops is very low density. 

We must accentuate the fact that, because for reaching an 
77* ~ 1 we will need an magnetic field of around 300 gauss for 
( Pup, a typical O-star with 5OM , while for the solar wind we 
have an 77* ~ 40. So, not just because of the different accel- 
eration mechanism, we can not account for the same topology 
of the magnetic field and for its effects on the wind particles 
for the Sun or for an O-type star, especially on the same phase 
space distribution. 

It must be emphasized that the Luminous Blue Variables 
(LBVs) giant outbursts may occur when these massive stars 
approach their Eddington limits. When this happens, they must 
reach a point where the centrifugal force and the radiative ac- 
celeration cancel out gravity at equator (Q-limit). When stars 
are close to the Q-limit, highly non-spherical mass loss should 
occur. This suggests a scenario where a slow and very dense 
wind, strongly confined to the equatorial pla ne, is fo llowed 
by a fast and almost spherical wind (Lang eretal. 19991) . With 
the assumption that the star is very close to the O-limit Langer 
et al. dLanger et al. 19991) models produced gas distributions 
that strongly resemble Homunculus, the hourglass-shaped in- 
ner part of the highly structured circumstellar nebula. 



CPup, ri.=3 2, lowMdot iso 

s "Ok sec 




Fig. 3. Final state for f Pup (M = 50M o ; R = 1.3 x 10 12 cm; 
M = 2.6 x 10~ 6 Mo/year; v M = 2500 km/ s) with rotating wind, 
77, = 32 and oj = D/Q c = 1 /2simulation Jud-Doula 20041) . The 
isothermal models with VTO < 77* < 10 and no rotation show a 
similar final state behavior 

The above models will consist our basic assumption for 
pre-supernova ejecta as a bipolar, dense, equatorial outflow. 

5.2. Shock Acceleration and Spallation of Cosmic 
Rays 

As was proposed by Biermann (Biermann 1993h and in 
Rachen & Biermann llRachen & Biermann 19931) . the cosmic 
rays observed events can be traced back to three sources 
llWiebel-Sooth et al. 19981) : 

1 . Supernova explosions into the interstellar medium, 
or ISM-SNe. This component produces mostly hy- 
drogen and the observed energetic electrons up to 
about 30 GeV, and dominates the all particle flux 
up to about 10 4 GeV. 

2. Supernova explosion into predecessor stellar wind, 
or wind-SNe. This component produces the ob- 
served electrons above 30 GeV, helium and most 
heavier elements already from GeV particle ener- 
gies. Due to a reduction in acceleration efficiency 
at a particular rigidity the slope of the spectrum in- 
creases, thus producing the knee feature. The com- 
ponent extends ultimately to several EeV Since the 
winds of massive stars are enriched late in their life, 
this component shows a heavy element abundance 
which is strongly increased over that of the ISM. 

3. Powerful radio galaxies produce a contribution 
which dominates beyond about 3 EeV, and consist 
mostly of hydrogen and helium, with only little ad- 
dition of heavy elements. At this energy the inter- 
action with the microwave background cuts off the 
contribution from distant extragalactic sources, the 
Greisen-Zatsepin-Kuzmin (GZK) cutoff. There are 
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a small number of events which appear to be beyond 
this energy, and whether they fit into such a picture 
is open at present. 

Various tests for these models were performed i n 
Biermann & Ca ssinelli dBier mann & Cassinelli"T9 93). 
Biermann & Strom (lBiermann"&^trorn*T993lT Stanev et al. 
( Stan ev etal. 19931) etc. 

The Cosmic Rays spectrum at GeV energies is close to 
E~ 2n5 , and at higher energies close to E~ 2 65 below the "knee" 
at » 5 • 10 15 eV, from where the spectrum turns down to about 
E~ 31 , to flatten out again near 3 • 10 18 eV (the "ankle"). The 
chemical composition is roughly similar to that of the inter- 
stellar medium, with reduced hydrogen and helium relative to 
silicon, and with some general enhancement of elements of low 
first ionization potential as we find in solar energetic particles. 

The theory which claims to account for the origin of cosmic 
rays below 3 • 10 18 eV traces them to the shock waves caused 
by supernovae exploding either into the interstellar medium 
(cosmic rays with energies bellow the "knee"), either into the 
predecessor stellar wind (cosmi c rays with energies above the 
"knee" but bellow the "ankle") dPopescu et al. 20051) . 

Into the picture of a nonspherical wind as above, we will 
expect that, when the supernova shock will hit the material sur- 
rounding the star it will accelerate particles to cosmic ray ener- 
gies through first order Fermi acceleration mechanism. At the 
magnetic equator of our massive star, due to the dense wind 
presence, it will be possible to accelerate particles to ~ 10 18 
eV. This acceleration limit will decrease dramatically from the 
equator to the polar regions because there practically the shock 
propagates in the ISM as it moves further away from the stellar 
surface. So, due to the time spent in shock (considered spher- 
ical), the maximal energy, which we will denote E, to which 
the particle will reach will be ~ 10 18 only in a narrow region 
around the equatorial plane, decreasing to around ~ 10 15 in 
the polar region. At low energies the isoenergetical surfaces for 
the shock accelerated particles will follow closely the isoden- 
sity surfaces (see Fig. This situation changes for energies 
10 15 < E < 10 18 eV, where the isoenergetical surfaces are fol- 
lowing the isodensity profiles to around < 30° in stellar latitude 
after which are cut off. The resulting geometry are catenoid- 
like surfaces centered on the star and having the stellar mag- 
netic axis as the rotation symmetry axis. 

The spallation is one of the processes that has one 
of the most important roles in the modification of 
the abundance of cosmic ray s in their trans portation 
through ISM dReeves et al. 1970t ISilberberg & Tsao 1973al 
Silb erberg & Tsao 1973bl) . The basic features of target frag- 
mentation, sometimes called "spallation", are very well 
understood: heavy fragments arise from peripheral collisions 
of heavy ions or relativistic protons wi th the target nucleu s 
dTsao et al. 19981 ISilberberg et al. 19981 ISihver et al. 19931 . 
These so called "spectators" of the reaction are excited 
primary fragments which then decay into the final fragments 
by a sequence of evaporation steps. 

The spectrum of the residual nuclei seems to be determined 
to large extent, but not fully, by the evaporation process. 



The difference between the spallation in high-energy 
physics in accelerators and in astrophysics is that, in the first 
case, the targets are the nuclei and the projectiles are the high- 
energy protons and in the second case the targets are the ISM 
protons (in first approximation) and the projectiles are the nu- 
clei that form the cosmic radiation. So everything depends just 
on which particle is put into the reference frame. 

About 10% by number of the interstellar gas is helium. 
Hence, about 20% of cosmic ray generated nuclear spal- 
lation nuclear products are formed in nucleus-helium inter- 
actions. In the analysis of cosmic ray interactions with at- 
mospheric nuclei (mainly nitrogen and oxygen) an accurate 
knowled ge of nu cleus-nucleus interactions becomes essential 
JTsao et al .19981). 

High-energy protons cause many different nuclear reac- 
tions and, in principle, all these nuclear processes have to be 
taken into account for a proper description of the total process. 
Possible outgoing particles are, for example, all light particles, 
gamma rays, neutrinos, and (above incident energies of 150 
MeV) pions. Furthermore, high-energy fission may occur and 
outgoing particles from the proton bombardment stage further 
reactions. Clearly, a wide spectrum of reaction products will be 
formed. 

The spallation happens when the supernova shock travels 
through the stellar wind and then hits the surrounding molec- 
ular shell of dense gas being responsible for the production of 
Li, Be and B in supernova shell. 

As well, the spallation with the ISM's protons or with the 
supernova predecessor star wind particles generates as sec- 
ondary particles pions (but also kaons and muons) which de- 
cay afterwards into muons and neutrinos secondary particles 
dCaisser 1990h: 

ff* -> if + v M (v p ) 
/f -> e ± + v e (v e ) + v^(v^) , 
with a similar chain for charged kaons. Also: 

p + A -> 7T° + B 



6. Concluding Remarks 

An essential feature of the gravitational interaction is that at 
center-of-mass energies well above the fundamental scale, the 
gravitational coupling grows so large that graviton exchange 
dominates over all other interactions. 

As said befor e, cosmic rays can probe the existence and 
the scale dAlvare z-Muniz et al. 2002[ |Fen^ & Shap ere 2002t 
lAnchordoaui & Goldberg 2002HEmparan et al. 20021) of extra- 
dimensions by observing the showers generated in Earth's 
atmosphere by brane decay. Black Holes decay ther- 
mally, according to the number of degrees of free- 
dom available, and so their decays are mainly hadroni c 
dGiddings & Thomas 2002t iDimopoulos & Lan dsberg 2 00 ll) . 
For example, the non-perturbative Black Hole formation and 
decay by Hawking evaporation gives rise to hard jets and lep- 
tons, with a characteristic ratio of hadronic to leptonic activity 
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Fig. 4. Catenoid 

of roughly 5: 1 . The high energy Black Hole cross section grows 
with energy at a rate determined by the dim ensionality and ge- 
ometry of the extra-dimensions (Giddinas & Thomas 2002). 

In this way, moving these physical arguments from Earth's 
atmosphere in the expanding supernova shell, even that the sec- 
ondary neutrinos produced in spallation are in MeV energy 
range, if the fundamental scale is 1 TeV, their scattering on 
the 10 18 eV accelerated particles may give rise to gravitons 
( V* ~ 10 5 GeV). Because the isoenergetical surfaces for these 
energies are catenoid-like surfaces, also the branes will form 
on a catenoid-like shaped surface, around the stellar remnant. 

Acknowledgements. I want to express all my gratitude to Valeriu 
Tudose whose comments regarding this paper's clarity where very 
useful to me and, I belive, also to the reader. 
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